Determining neutrino absorption spectra at Ultra-High Energies. 
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A very efficient method to measure the flux of Uitra-high energy (UHE) neutrinos is through 
the detection of radio waves which are emitted by the particie shower in the iunar regoiith. The 
highest acceptance is reached for radio waves in the frequency band of 100-200 MHz which can be 
measured with modern radio telescopes. In this work we investigate the sensitivity of this detection 
method to structures in the UHE neutrino spectrum caused by their absorption on the low-energy 
relic anti-neutrino background through the Z-boson resonance. The position of the absorption peak 
is sensitive to the neutrino mass and the redshift of the source. A new generation of low-frequency 
digital radio telescopes will provide excellent detection capabilities for measuring these radio pulses, 
thus making our consideration here very timely. 



I. INTRODUCTION 

With the advent of a new generation of digital ra- 
dio telescopes working in the frequency range of 100- 
200 MHz it will become possible to measure the flux of 
Ultra-High Energy (UHE) neutrinos, energies in excess of 
10 20 eV, as they impact the Moon PJ. To determine this 
flux as function of their energy is of considerable interest 
as it may give some clues about their origin. The recent 
results of the Pierre Auger Observatory on the correla- 
tion of UHE cosmic rays with the positions in the sky of 
AGN have only increased the interest. 

UHE neutrinos will be absorbed on the low energy Cos- 
mic Neutrino Background (Ci/B) through the excitation 
of the Z-boson resonance [3|. In general this annihila- 
tion process is expected to lead to sizeable absorption 
dips in the UHE neutrino spectra [H, 0] ■ The positions 
of these absorption lines depend on the redshift of the 
source emitting the UHE neutrinos. The issue addressed 
here is the precision with which these absorption dips 
can be determined in a realistic experiment. For the 
UHE background we will assume the flux as given by 
the Waxman-Bahcall (WB) limit Q. 

An exciting possibility is that the neutrino mass is not 
constant but varies with redshift, Mass Varying Neutri- 
nos (MaVaN). This has been proposed Q as a model 
for dark energy. Such a neutrino mass variation intro- 
duces a distinctly different structure of the absorption 
lines due to resonant absorption on relic anti-neutrinos 
into Z-bosons than is obtained for redshift independent 
masses 0, 0| • I n this work we will therefore also inves- 
tigate if this difference is observable in a realistic exper- 
iment. A very promising method [l[ for measuring the 
flux of UHE neutrinos is through the detection of the ra- 
dio pulse which is emitted when UHE neutrinos initiate 
a shower in the crust of the Moon. 

In Section |TT] a short quantitative discussion of radio 
emission from neutrinos hitting the lunar surface is pre- 
sented. Following the approach of Ref. [l[ we will calcu- 



late the pulse-height spectrum, as may be measured at 
Earth. First we investigate in Section fill Al the detector 
response for mono-energetic neutrinos. In Section IIII Bl 
the response to a spectrum of neutrino energies is given 
using different scenarios for the neutrino mass and red- 
shift of the source. Using the statistics of a possible ob- 
servation we show to what extent the spectra for different 
neutrino-absorption scenarios may be distinguished using 
a new generation of low- frequency digital radio telescopes 
such as LOFAR Q and SKA [jj]. The conclusions are 
presented in Section [TV] 
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II. RADIO EMISSION 

When an UHE particle induces a shower in dense me- 
dia, such as ice, salt, and lunar regoiith, the front end of 
the shower has a surplus of electrons. Since this cloud of 
negative charge is moving with a velocity which exceeds 
the velocity of light in the medium, Cerenkov radiation 
is emitted. For a wavelength of the same order of magni- 
tude as the typical size of this cloud, which is in the radio- 
frequency range, coherence builds up and the intensity of 
the emitted radiation reaches a maximum. This process 
is known as the Askaryan effect [l(| ■ Dagkesamanskii and 
Zheleznyk [Tlj were the first to apply this mechanism to 
UHE neutrinos hitting the moon and proposed to use ra- 
dio telescopes at Earth to detect the emitted radiation 
which is coherent for radio waves. Several experiments 
have since been performed [12|, |_13f| to find evidence for 
UHE neutrinos. All of these experiments have looked for 
this coherent radiation near the frequency where the in- 
tensity of the emitted radio waves is expected to reach 
its maximum. Since the typical lateral size of a shower is 
of the order of 10 cm the peak frequency is of the order 
of 3 GHz. Recently it has been proposed [l[ to perform 
such observations at an even larger wavelength of about 
3 m. This corresponds to the longitudinal extend of the 
shower in the lunar regoiith and as a result the angular 
spread of the emitted signal is close to maximal while 
keeping coherence. As a result, the detection efficiency 
for UHE particles is many orders of magnitude higher at 
100 MHz than at 3 GHz. 
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The intensity of radio emission (expressed in units 
of Jansky's where 1 Jy = 10 -26 W m _2 Hz _1 ) from a 
hadronic shower, with energy E s , in the lunar regolith, 
in a bandwidth Ais at a frequency v and an angle can 
be parameterized as [l| 



F(6,v,E s ) = 3.86 x 10 4 / 



(d>moon \ ( 



sin6» c y VlO 20 eV 

2/ A!/ 



H>(1 + (^o) 1 - 44 ) 



100 MHz 



with 
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Vn 2 - 1/ VttA c 



Jy,(i) 
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where i>o = 2.5 GHz [IsT ]. d is the distance to the observer, 
and d moon = 3.844 x 10 8 m is the average Earth-Moon 
distance. The angle at which the intensity of the radia- 
tion reaches a maximum, the Cerenkov angle, is related to 
the index of refraction (n) of the medium, cos# c = l/n. 
The parametrization of the energy dependence of the 
spreading of the radiated intensity around the Cerenkov 
angle, A c , is given in Rcf. [![. 

The mean value for the attenuation length of radio 
waves in the regolith is taken as X r = (9/i/[GHz]) m. In 
the calculations we have included radiation coming from 
a depth of at most 500 m treating for simplicity, and 
without loss of accuracy [l|, the whole layer as behaving 
like regolith. For the neutrino an energy-dependent mean 

1/3 

free path [3 has been used, \ v = 130 ( 10 ^ cV ^l km, 



which is appropriate for regolith. 

For neutrino-induced showers only 20% of the initial 
energy is converted to a hadronic shower, while the re- 
maining 80% is carried off by the lepton. This ener- 
getic lepton will not induce a detectable radio shower. 
For a muon, the density of charged particles will be too 
small, while the shower of an UHE electron will be elon- 
gated due to the Landau-Pomeranchuk-Migdal (LPM) 
effect 0]. The width of the Cerenkov cone will thus 
be very small which makes the shower practically unde- 
tectable. For the present calculations we therefore have 
limited ourselves to the hadronic part of the shower. It 
should be realized that the photons in the hadronic part 
of the shower have an energy that is roughly 2 orders of 
magnitude less than the primary energy and thus that 
the LPM effect is considerably smaller for this part. 

The radio flashes from the moon can be measured with 
powerful synthesis telescopes such as LOFAR [§]. LO- 
FAR can operate in the frequency band of 115-240 MHz 
where it will have a sensitivity of about F noisc =20 Jy. 
Simulations show that a pulse of intensity 25xF' noiso , in- 
terfering with the noise background, can be detected with 
3a significance at a probability greater than 80%. For 
this reason we have assumed in the calculations a detec- 
tion threshold of 25xF noiao . A simulation for LOFAR, 
taking v = 165 MHz, a bandwidth of Ai> = 100 MHz, 
and an observation time of 30 days is shown in Fig. [1] 
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FIG. 1: Flux limits, summed over all flavors, on UHE neutri- 
nos as can be determined with LOFAR observations (see text) 
are compared with various models, in particular, WB [j (ver- 
tical bars) , GZK [TiT] (dotted thin green line). Limits from the 
RICE pi, ANITAliai, and FORTE [ij experiments are also 
shown. 



assuming a 60% Moon coverage. The result is compared 
with the WB limit Q . With the future SKA telescope a 
sensitivity can be reached 0] which is higher by several 
orders of magnitude. 

The following analysis of the spectral sensitivity of the 
method is in principle independent of the threshold sensi- 
tivity for detecting a radio signal. However with a higher 
sensitivity a larger fraction of the pulse-height spectrum 
can be determined with the possibility of increasing the 
spectrum-resolving power. 



III. SPECTRAL RESOLVING POWER 

An appealing opportunity to catch a glimpse of the 
CuB is through the energy spectrum of UHE neutri- 
nos. Such UHE neutrinos can annihilate with relic anti- 
neutrinos (and vice versa) into Z bosons, if their energies 
coincide with the respective resonance energies E^f of 
the corresponding process v + v — > Z . These energies, 



77 res 



2m M 



= 4.2 x 10 



21 



eV 



-eV 



(3) 



in the rest system of the C^B, are entirely determined 
by the Z boson mass Mz as well as the respective neu- 
trino masses m„ i . An exceptional loss of transparency 
of the UHE neutrinos results from the fact that the cor- 
responding annihilation cross-section on resonance is en- 
hanced by several orders of magnitude with respect to 
non-resonant scattering. As a consequence, the diffuse 
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UHE neutrino flux arriving at Earth should show absorp- 
tion lines due to this resonance scattering off the CVB, as 
explained in detail in Ref. [3]. Our aim here is to deter- 
mine the precision with which it is possible to measure 
such absorption dips. As reference or source spectrum, on 
top of which the absorption is imposed, the WB limit Q 
is taken, 

<$>wb = E 2< ^ = 3.5 x 1CT 8 [GeV cm^s^sr- 1 ] . (4) 
dE 

It is not a straightforward issue to determine the original 
neutrino energy spectrum from the measured signals at 
Earth. When a neutrino interacts in the lunar regolith, 
the magnitude of the generated electric pulse, as it arrives 
at Earth, is only indirectly related to the original energy, 
as can be seen from Eq. JT]). The strength of the pulse 
depends strongly on emission angle and the energy frac- 
tion of the hadronic part of the shower [2l|. In addition 
the signal is attenuated in the lunar regolith. Although 
it may be possible to determine emission angle and depth 
from the frequency dependence of the pulse, we do not 
want to rely on this in the present investigation. To cal- 
culate the sensitivity of the NuMoon measurements to 
the energy spectrum of UHE neutrinos, as a first step the 
pulse-height spectra are calculated as received at Earth 
from mono-energetic neutrinos, integrated over angles of 
incidence, emission angles, and emission depth. In the 
second step these pulse-height spectra are integrated over 
the neutrino-energy distributions. 

A. Single energy 

To be able to address quantitatively the pulse-height 
spectrum which is measurable at Earth, we introduce 
the detection probability T>(h;E). It is defined as the 
probability that a neutrino hitting the Moon at an arbi- 
trary angle and position with energy E will produce radio 
waves at a frequency of v = 100 MHz, with a power at 
Earth exceeding a certain threshold value ft, in a band- 
width of 50 MHz. In Fig. [2] detection probabilities are 
given as a function of the detection threshold h for dif- 
ferent neutrino energies in the range of interest for this 
paper. 

The probability V(h; E) for measuring a signal with a 
pulse height between ft, — 0.5 and ft + 0.5 when a neutrino 
with energy E impinges on the surface of the Moon, can 
now be obtained by differentiating T>(h; E) with respect 
to the detection threshold h. The obtained functional 
dependence of V(h; E) on ft is shown in Fig. [3] 

The maximum pulse height which can be detected at 
Earth emitted by a moon shower which is induced by a 
neutrino with an energy E is written as ft = A(E). This 
maximum follows from the strength of the emitted radio 
pulse as given by Eq. (P) plus some 20% attenuation at 
the lunar surface, 

A(E) = 128 Jyx (j^v) , (5) 



T 




h[Jy] 

FIG. 2: Detection probability T>(h; E) v.s. threshold h for dif- 
ferent neutrino energies E. Blue, dotted, line: E = 10 21 eV, 
magenta, dashed, line: E — 2 x 10 21 eV, black, close-dotted, 
line: E = 4 x 10 21 eV, red, drawn, line: E = 10 22 eV, brown, 
close-dashed, line: E = 10 23 eV. 

and this reproduces the end-point values given in Fig. [5] 
Of particular interest here is the energy dependence 
which will be useful to parameterize the energy de- 
pendence of P(h;E). Numerical results indicate that 
P(h; E) scales with energy as 

(E'\ 1 ' 73 
~e) (6) 

with 

A(E) \E J v ' 

where A(E) is given by Eq. (J5J). To show the accuracy of 
this parametrization we plot in Fig. [4] the ratio of the ex- 
act probabilities at E=10 22 eV with those reconstructed 
from the results at lower energies. Considering the sim- 
plicity of the parametrization over several orders of mag- 
nitude, the agreement is excellent. 

B. The neutrino absorption spectrum 

The absorption of UHE neutrinos in the intergalac- 
tic space is generally small and the neutrinos arriving 
at Earth may have originated from sources at great dis- 
tances. In the calculation of the absorption dips in the 
spectrum it is therefore essential to include in the consid- 
eration the expansion of the universe and thus the red- 
shift of far-away sources. In the following we discuss the 
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FIG. 3: Probability distribution P(h;E) for neutrinos with 
the same energies as used in Fig. [2] 



structure of the absorption dips in the neutrino spectrum 
in two different pictures, first for the picture in which the 
neutrino mass is independent of redshift followed by that 
for the MaVaN picture Q. 

For a redshift independent neutrino mass the width of 
the absorption peaks depends strongly on the redshift of 
the source, z s . The peak position is determined by the 
neutrino mass. Without considering thermal broadening, 
which is a small effect, resonant absorption occurs for the 
neutrinos in the energy interval Q 

-Ji- < E < Eft (8) 
1 + z a 

where E^ e f is given in Eq. ([3]). Absorption at energies be- 
low E^f occurs at large distances where the neutrino en- 
ergy is larger than observed at Earth due to the redshift. 
The thick, blue, curves in Fig.[5]give the neutrino absorp- 
tion spectra, using the values given in Fig. 5 of Ref. [fj], 
in the picture that the neutrino masses are independent 
of redshift, the static-mass neutrino (StaMaN) picture. 
The results are given for three different distances of the 
source, corresponding to redshifts of z s = 5, z s = 20 and 
z s = 50. A normal hierarchy for neutrino masses has 
been used, m u i = 10~ 5 eV, m V 2 = 8.3 x 10~ 3 eV and 
mv3 = 5.17 x IO -2 eV. The absorbtion lines are calcu- 
lated for a universe with a day matter density Qm = 0.3, 
a curvature density flk = and vacuum energy den- 
sity f^A = 0.7 as corresponds to a Lambda Cold Dark 
Matter (ACDM) universe as suggested by recent obser- 
vations [161 ] with a homogeneous relic neutrino density. 

The absorption lines are superimposed on a smooth 
neutrino spectrum with a flux equal to the WB-limit. 



FIG. 4: Ratio of the exact probabilities at E=10 eV with 
those reconstructed by applying Eq. [6] from the results at 
lower energies; blue, dotted, line from E'=10 21 eV, magenta, 
dashed, line from E'=2 x 10 21 eV and black, dashed-dotted, 
line from E'=4 x 10 21 eV. 



From Fig. [S] one sees clearly that for larger values of z s 
the absorption dips show at smaller neutrino energies, as 
predicted based on Eq. ©. In addition the absorption is 
getting more pronounced because of the larger travelling 
distance of the UHE neutrinos through the C^B. Struc- 
tures in Fig. [5] at lower energies shown are due to the 
absorption of neutrinos with higher mass. Due to neu- 
trino oscillations the absorption at lower energies can be 
almost complete, instead of only 1 /3 when only neutrinos 
of a single flavor are absorbed. Most important for the 
present study are the heavier mass neutrinos which gen- 
erate structures in the spectrum at lower energies where 
the sensitivity for detecting it is largest due to the higher 
flux. 

In the MaVaN model the absorption dips in the neu- 
trino spectrum become independent of z a for small z s 
because m ui (z) oc (1 + z)^ 1 . This implies that for rel- 
atively nearby sources of UHE neutrinos one expects to 
see sharp absorption lines in the spectrum, given by the 
thin, red, curves in Fig. [5] For sources at distances cor- 
responding to z s > 10 the predicted redshift dependence 
is m l/ i(z) oc (I + z) -1 / 2 and as a result the lines broaden 
considerably. At a redshift of z s = 50 this effect is so 
big that the differences in the absorption spectra in the 
MaVaN picture and the StaMaN picture are minor Q at 
smaller neutrino energies. 

Based on the detection probability for pulses at Earth, 
the neutrino flux, and the length of a typical experiment, 
the expected number of radio pulses of a particular mag- 
nitude can be calculated. For the present estimates given 
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FIG. 5: The averaged neutrino survival probability, 
1/3 X] Pvt , as a function of neutrino energy (extracted from 
Fig. 5 of Ref. [|j]) for sources at different redshifts; z a = 5 
(drawn line), z s = 20 (dotted line), z s = 50 (dashed line). 
The thick, blue, curves give the absorption spectra for the case 
that the neutrino mass is independent of redshift, while the 
thin, red, curves gives the predictions in the MaVaN model. 



in Fig. [5] the typical measuring conditions with LOFAR, 
as discussed in Ref. are assumed, i.e. a measuring 
time of 30 days where the telescopes have only half the 
moon is in the field of view. The quoted error bars reflect 
the statistical error only. Of course for longer observing 
times or better coverage of the moon the statistics will 
improve. The bin size of the pulse-height spectrum is 
chosen the same as was assumed in Fig. [3] where the size 
of the bin is taken equal to the pulse height . 

From Fig. [6] it can be seen that in a realistic experi- 
ment one is sensitive to the absorption structures in the 
spectra only if the sources are at the largest redshifts. 
Even though the absorption structures in the neutrino 
spectra are distinctly different, these structures are close 
to impossible to unfold from the measured pulse-height 
spectrum. One thus concludes that, on the basis of these 
data, it will not be possible to distinguish the MaVaN 
and the StaMaN pictures. On the other hand the mea- 
sured values are clearly sensitive to the source distance, 
independent of the picture used for the neutrino mass. 

For a realistic observation, however, the absolute neu- 
trino flux in not known and the absolute number of 
counts can not be used to distinguish the different fea- 
tures in the neutrino spectra. For this reason we com- 
pare the spectra where the counts have been normal- 
ized to the same value at a relatively low pulse height, 
h = 4 x F aoiso — 80 Jy in the present example, which 
should be typical for a LOFAR observation. The normal- 
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FIG. 6: The expected count rate for a realistic experiment for 
the different neutrino absorption dips given in Fig. [5] 
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FIG. 7: The counts given in Fig.|6]are re-normalized to those 
at a pulse height of 80 Jy to the value expected from the 
spectrum without absorption. 

ized spectra are shown in Fig. [7l From this figure it can 
be seen that the observation will not be able to convinc- 
ingly distinguish between the MaVaN and the StaMaN 
models. One should however be able to determine the ab- 
sorption lines if the sources are at distances of the order 
of z s = 50. Seeing the absorption lines would confirm 
that indeed a Ci^B exists. In observations with SKA, 
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where the higher sensitivity allows to distinguish even 
smaller pulses from the background, one may be able to 
see absorption structures with sources as close as z s = 20. 

IV. SUMMARY 

Due to resonant absorption on the CjaB the energy 
spectrum of UHE neutrinos carries with it the informa- 
tion on the neutrino mass and the distance of the emit- 
ting source from the Earth. As stressed already in Ref. [f| 
the structure of the absorption lines is markedly differ- 
ent in a model where the neutrino mass is constant as 
compared to that predicted in the MaVaN model @. In 
this work we have addressed the question whether these 
differences can be determined in a realistic experiment. 
The most realistic possibility to observe the small flux 
expected for UHE neutrinos is by observing the radio 
waves in the frequency range of 100-200 MHz that are 
emitted when these neutrinos impact the Moon pj . Since 
the observed pulse height is only indirectly related to the 



primary energy of the neutrino we have made a realistic 
simulation of such a measurement based on the neutrino 
absorption calculated in Ref. 0. The results show that, 
due to the smallness of the absorption or, equivalently, 
due to the high degree of transparency of intergalactic 
space for UHE neutrinos, absorption lines arc hardly vis- 
ible for sources closer than z s = 20. For sources at a 
redshift of z s = 50, the effects of the absorption on the 
CVB is clearly visible in the measurement. The predicted 
differences between the two models for the neutrino mass 
are too small to be detectable in a convincing way. 
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